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1. Introduction 

Within last century the study of cycles of two- 
dimensional quadratic systems was stimulated by 
16th Hilbert's problem and its different variants 
[Hilbert, 1902; Lloyd, 1988; Blows & Perko, 1994; 
Ilyashenko, 2002; Chen & Wang, 1979; Shi, 1980]. 

The development of asymptotical methods for 
investigation of cycles of quadratic systems is fa- 
cilitated by the fractionally-quadratic transforma- 
tions, which reduce any two-dimensional quadratic 
system to the Lienard equation of special type 
with discontinuous nonlinear functions [Leonov, 
1997; Leonov, 1998; Leonov, 2006]. In the present 
paper these new facilities are demonstrated. 

The elements of direct Lyapunov method, the 
small perturbations methods, and the bifurcation 
theory together permits us to obtain the existence 
criteria of small cycles in the neighborhoods of sta- 
tionary points of Lienard equation. 

This criteria applied to quadratic systems. 

The question arises as to whether in two-dimen- 
sional quadratic systems a simultaneous bifurca- 
tion of cycles from two weak stable and weak un- 
stable equilibria is possible by variation of one 
scalar parameter. Here we provide the positive 
answer to this question. 

For the Lienard equation, the classical theorems 
on the existence of cycles is well known [Lefschetz, 
1957; Cesari, 1959]. Here there occurs quite nat- 
urally the possibility of generalization of these re- 
sults and of applying them to the problem on the 
existence of cycles in quadratic systems. On this 
way in the space of parameters of quadratic sys- 
tems, the sets of positive Lebesgue measure, for 
which the cycles exist, are discovered. 

2. The transformation of quadratic system 
to the Lienard equation 



Let us consider the quadratic system 
x = a\x 2 + b\xy + c\y 2 + a±x + P\y 
y = a 2 x 2 + b 2 xy + c 2 y 2 + a 2 x + /3 2 y 

where dj, hi, Ci,oti, Pi are real numbers. 

Here we shall follow the paper [Leonov, 1997; 
Leonov, 1998; Leonov, 2006]. 

Proposition 1. Without loss of generality, we 
can assume that c\ = 0. 

Proof. Suppose, for definiteness, that a 2 ^ 
(otherwise, having performed the change x — > y, 
y — > x, we obtain at once c\ = 0). Introduce 
the linear transformation x\ = x + uy, y\ = y. To 
prove Proposition 1, it is sufficient to show that for 
the certain numbers p, k, v the following identity 

(x -I- vyf = (ai + ua 2 )x 2 + (b\ + ub 2 )xy+ 

+ (ci + vc 2 )y 2 + (ai + va 2 )x + (Pi + v/3 2 )y = 

= p(x + vy)y + k(x + vy) 2 + {ct\ + va 2 )x+ 

+ (Pi + vfa)y 

is valid. This identity is equivalent to the following 
system of equations 

k = a\ + ua 2 , 

kv 2 + pv = ci + uc 2 , (2) 
p + 2nu = b\ + vb 2 . 

The above relations are satisfied if 

(ai + va 2 )v 2 - v(b\ + vb 2 ) + (ci + uc 2 ) = 0. 

Since a 2 ^ 0, this equation with respect to u al- 
ways has a real root. Thus, the system of equa- 
tions (2) always has a real solution. □ 
Further we assume that c% = 0. 

Proposition 2. Let be b\ ^ 0. The straight 
line Pi + b±x = on the plane {x, y} is invariant 
or transversal for trajectories of system (1). 



Proof. The last statement follows from the re- 
lation 

(Pi + b 1 x)'=b 1 [(bix + (3i)y + a x x 2 + a\x\ = 
2 



1 /T 



.ft 



for x = —Pi/bi. This implies that if 

then the straight line Pi + b\x = is invariant and 
if 



ai 



Pi 



n , | ^ ] - 0. 



then the straight line /3i + &ix = is transversal. 
□ 

Excluding the trivial case that the right-hand 
side of the first equation of system (1) is indepen- 
dent of y, we assume that 

l&il + lftl^O. (3) 

Then by Proposition 2 we obtain that the tra- 
jectories of system (1) are also the trajectories of 
the system 



x = y + 



CL\X + U\X 

Pi + bix 



y 



a 2 x 2 + b 2 xy + c 2 y 2 + a 2 x + (3 2 y 



(4) 



Pi + hx 

Consider the following transformation 

ci\x 2 + a.\x 



y = y + 



x = x. 



Pi + b\x 



Using these new phase variables (here the bars 
over the variables x — > x, y — > y are omitted), 
system (4) can be written as 

x = y, y = -Q(x)y 2 - R(x)y - P(x), (5) 

where 



Q(x) = 



-c 2 



R(x) 



P(x) 



Pi + bix' 

(bib 2 - 2aic 2 + aibi)x 2 + (& 2 ffi + hp 2 - 
(Pi + M 2 

-2aic 2 + 2oi)9i)g + ai^i + Pip 2 
(Pi + M) 2 
a 2 x 2 + a 2 x (b 2 x + p 2 )(aix 2 + uix) 



Pi + bix 
c 2 (aix 2 +aix) 2 
(Pi+bixf 



[Pi + hx) 2 



Proposition 2 and condition (3) imply that the 
trajectories of system (5) are also trajectories of 
the system 

x = ye p(x \ y= [- Q(x)y 2 - R(x)y - P(x)]e pix \ 

where p(x) is a certain integral of the function 
Q(x). 

From this system, using the change x = x, y = 
we obtain 

x = y, y = -f(x)y - g(x). (6) 

Here the bars over the variables x and y are also 
omitted. 

So, for bi ^ 0, by means of the above nonde- 
generate changes we can reduce system (1) to the 
Lienard equation (6) with the functions 

f(x) = R(x)e p W = R(x)\pi + bix\ g , 

g(x) = P(x)e 2p{x ^ = P(x)\Pi + bix\ 2q 

Here for &i ^ 0, q = -g. 

Proposition 3. Let be bi ^ 0, Pi ^ 0, cci / 
0. Without loss of generality, we can assume that 
bi = ai = Pi = 1. 

Proof. Using the change 

Pi- oil- , t 
x = -7-x, y = —y, t = — , 
bi bi ai 



we obtain 



- aiPiai_ 2 

x = — x + xy + x + y. 



□ 



Let us consider system (6) where 
fix) = (Ax + B)x\x + l\ q - 2 , 
. (Cix 3 + C 2 x 2 + C 3 x + C A )x 2 

9(X) = j . TTq \X + 1 Q . 

(x + l) d 

Here A, B, Cj (j = 1, . . . , 4), q are real numbers. 

Proposition 4. For numbers A, B, Cj,q of sys- 
tem (6) exist corresponding numbers a±, &i = 1, cci 
1, Pi = 1, a 2 ,b 2 ,c 2 = -q, a 2 ,p 2 = -1 of system 
(l)iff 
(B-A) 



(2q - l) 2 
(B-A) 

(2g - 1; 

Here 



((1 - q)B + (3q - 2) A) = 2C 2 - 3d - C 3 , 
2 (B + 2(g - 1)A) = C 2 - 2Ci - C 4 . 

B-A 



ai = 1+ 2q-V 
0-2 = —(q + - ^4oi - Ci, 
b 2 = -A-ai(2q+l), 
a 2 = a\- 2ai + A(a x - 1) + (2Ci - C 2 ). 



3. Existence criterion for small cycles of 
the Lienard equation in the neighborhood 
of equilibrium 

Consider the equation 

z + z = u(t). (7) 

For this equation it is well known the following 
simple result. 

Lemma 1 [Arnol'd, 1976]. The solution of equa- 
tion (7) with initial data z(0),i(0) is given by for- 
mula 

t 



z(t) 



z(0) — J u(t) sin rdr 



+ 



cos rdr 



cost+ 



sint. 



(8) 



For i(0) = 0, taking into account formula (8), 
we have the following relation 



z(t) = — z(0) sin t + J u(t) sin rdr sin t+ 



+ 



J u(r) 



(9) 



cos rdr cost. 



Consider the equation 

x + F(x,e)x + G(x,e) =0 

or the equivalent system 
x = y 

y = -F(x,e)y - G(x,e). 



(10) 



Here e is a positive number, F(x, 0) = f(x), G(x, 0) 
g(x), f(xo) = g(xo) = 0, in a certain neighbor- 
hood of the point x = xq,£ = the functions 
F(x, e) and G(x, e) are smooth functions. 
Theorem 1. If the inequalities 

f"(x )g'(x )-g"(x )f(x )<0, 

g'(x o )>0, F(x £ ,e)>0 1 ' 

where x £ is a zero of the function G(x, e) in the 
neighborhood of the point x = xq, are satisfied, 
then there exists a number £o > such that for all 
e G (0, £o) system (10) has a cycle. 

Proof. We introduce the notation z(t) = x(t) — 
xo and put 

f(x) = h{x - x ) + f 2 (x - x ) 2 + 0((x - x ) 3 ), 
9{x) = gi(x - x ) + g 2 (x - xof + 0((x - x ) 3 ). 



Without loss of generality, it can be assumed that 
91 = 1- 

Consider the case when e = and z(0) is a small 
number. 

Here we shall use smoothnis of functions F and 
G and shall follow the first Lyapunov method on 
finite time interval [Lefschetz, 1957; Cesari, 1959]. 

The first approximation to the solution z(t) of 
system (10) is the function 

z\{t) = z(0) cost. 

The second approximation z 2 (t) to the solution 
z(t) is given by the equation 

z 2 + z 2 =u(t), (12) 

where 

u(t) = z(0) 2 (/i cost sin t — ^(cost) 2 ). 
By Lemma 1 we obtain that z 2 (t) has the form 

z 2 (t) = z(0)cos£ + z(0) 2 ^-y(sin£) 3 cos£+ 

+ ^(1 - (cost) 3 ) cost + ^(1 - (cost) 3 ) sin t- 

1 \ 

— g 2 (smt — -(sint) ) sint J . 

(13) 

This implies that for z(0) = the relation holds: 

z 2 (t) = -z(0)smt + z(0) 2 (J±(sint) 4 - 

— — (1 - (cost) 3 ) sint + —(1 - (cost) 3 ) cost- 

1 \ 

— g 2 (sin t — - (sin t) ) cos t j . 

It follows that the time T > 0, of the second 
crossing, of the straight line {y = 0} by the solu- 
tion x(t),y(t) of system (10) with the initial data 
x(0), y(0) = satisfies the relation 

T = 2^ + O((x(0)-x ) 2 ). (14) 
Consider now the function 

/ x \ 2 x 

V(x,y)= ly + J f(z)dz) +2 J g(z)dz. 



X'() 



For the derivative of the function V along the 
solution of system (10) the following relation 



V(x,y) = -2g(x) J f(z)dz 



= -h{x - x ) 3 - l-f 2 + f ig2 ) (x - x ) 4 + 

+ O((x-x ) 5 ) 

is valid. Then, taking into account relations (13), 
(14), we obtain 

V(x(T),y(T))-V(x(0),0) = 

T 

= ~ j (fl^iff + f |/ 2 + /lift) Z 2 (t)^J dt+ 



+ O(z(0) 5 ) = -2(0)' 



2vr 

J {^h + hg^j (cost) 4 + 



+ 3/i(cost) 2/ ^ Lf - 



92, 



(sint) d cost + ^(1- 

t 

(cost) 3 ) cost + —(1 - (cost) 3 ) sin t- 
3 

i \ 

52 (sin t— — (sint) )sintl dt+ 



+ O(z(0) 5 ) 



z(0) A + O(z(0f). 



In this case by the theorem on a continuous depen- 
dence of solutions of system (10) on a parameter, 
for sufficiently small parameter e > in compari- 
son with the number \z(0)\ we have the inequality 

V(x(T),y(T))>V(x(0),0). 

On the other hand, for small e > the equilib- 
rium x = x E , y = of system (10) is a stable focus. 
These two facts yield that in a certain small neigh- 
borhood of the point x = x e ,y = there exists a 
cycle of system (10) (Fig. 1). Theorem is proved. 
□ 

Example 1. Consider the equation 

x + (e + fix + f2X 2 )x + x + g 2 x 2 + # 3 x 3 = 0, 

where fx, f 2 , g 2 ,g 3 are arbitrary numbers, / 2 -/i52 < 
0. 

Theorem 1 implies that for the sufficiently small 
positive £ this equation has at least one periodic 
solution. 




Figure 1. 

Consider system (1) and assume that 

h = fh = 1, «i + fa = e, 

A = -b 2 + 2a\c 2 - ai, 

B = —b 2 —02 + 2aic 2 - 2oi, 

C = —a 2 - 2a 2 + a\b 2 + a\f3 2 + otxfc ~ c 2 a\, 

D = —a 2 + ai/3 2 - 

If £> > and 

AD - BC + BD(1 + c 2 ) <0 

then by theorem 1 we obtain that for sufficiently 
small e > the system (1) has a cycle in some 
neighborhood of the point x = y = 0. 

Suppose that in system (1) a\ = a 2 = c\ = 0, 
h = Pi = c 2 = 1, 6 2 = -1, ai = 1/3 — e, a 2 = 
(3 2 = -1/3. We have 

F(x,e) = (x 2 + 2(l-e)x + e)|l + x|~ 3 , 

G(x,e) = [ \{x 2 + 2x) + ex(x 2 + 2x + - - £x] 
\ 9 3 

(1 + x)- 5 , 



/'(0) = /'(-2) = 2, </(0) = </(-2) = - 
f(0) = /"(-2) = -10, /(0) = /(-2) 



-2. 



In the neighborhood of x = 0, x £ = and in 
the neighborhood of x = —2 the relation x £ = 
-2 - 3£ + o(e) holds. 

For x e = 0, we have 

F(x £ ,e) = £ 

and for x e = — 2 — 3e + o(e), 

N ll£ + o(e) 
1 + x e r 



It should also be noted that in the considered 
case system (10) with e = under the change 

x = —z — 2, y = —i? 

preserves its form and the equilibrium x = —2, y = 
goes into the equilibrium z = •& = 0. 

Thus, we prove the following 

Theorem 2. For sufficiently small e < 0, the 
system 

1 



x = xy+[--e)x + y 

2 1 1 

y = -xy + y - -x - -y 

has no less than two cycles. Moreover, each of two 
equilibria is surrounded by no less than one cycle. 



4. The generalization of the Lienard 
theorem 

In this section we give the generalization [Leonov, 
2006] of a classical result of Lienard concerning 
the existence of periodic solution for the equation 
[Lefschetz, 1957] 

x + f(x)x + g(x) = (15) 

Consider a system 

x = y 

v = -f{x)y-g{x)' 



(16) 



which is equivalent to equation (15). 

Suppose, the function f(x) and g(x) are contin- 
uous on the interval (a, +oo) and for the certain 
numbers a < v\ < xq < v 2 the following condi- 
tions 

1) lim g(x) = — oo, 

lim g(x) = +oo, 

x^+oo 

X X 

lim / g(z)dz = lim / g(z)dz = +oo, 

x^a J x^+oo J 



Xo 



Xo 



2) f(x)>0, VxG (a,^i)(J(i/ 2 ,+oo), 



J f(x)dx 



< 



(17) 



(18) 



V2 

are satisfied. 

Theorem 3. // conditions 1) and 2) are valid, 
then in the phase space 



there exists a piecewise- smooth transversal closed 
curve, which intersects the straight line {y = 0} at 
the certain points a < n\ < v\ and fj, 2 > ^2 (Fig. 
2). If, in addition, in space (19) system (16) has 
only one unstable by Lyapunov focal equilibrium, 
then system (16) has a cycle. 

Proof. Consider a pair of the numbers fii G 
(a,vi) and fi 2 G (u 2 , +00) such that fi\ is suffi- 
ciently close to a, fi 2 is sufficiently large, and 



f*2 



g{x)dx = 0. 



(20) 



mi 



Without loss of generality, we put 
g(x) < 0, \/x G [ni,v{\, 

g{x) > 0, Mx G [v 2 ,fJ, 2 ]. 
Consider the following functions 



V 1 (x,y)=y 2 + 2 J g(z)dz, 

xo 

(x \ 2 x 

y + j f(z)dz\ +2 J g(z)dz, 
vi I xo 

/ x \ 2 x 

V 3 (x,y) = ly + Jf(z)dz\ +2jg(z)dz, 



V2 



Xo 



{x G (a, +00), y G -R 1 } 



(19) 



V A (x,y) = V 2 (x,y) - e(x - vi) 
V${x,y) = V 3 (x,y) - e(x - v 2 ) 
Ve{x,y) = V 3 (x,y) + e(x - v 2 ) 
V7{x,y) = V 2 (x,y) + e(x - v\). 

Here e is a certain sufficiently small number. 
Consider now the sets Qj (see Fig. 2) 



[fj,i,vi],y > 0,Vi(x,y) < ¥l(^i,0)}, 
[vi,x ],y > Q,V 4 (x,y) < V 2 (v u yi)}, 
[x ,u 2 },y > 0,V 5 (x,y) < V 3 (u 2 ,y 2 )}, 
[v 2 ,n 2 ],y > 0,V 3 (x,y) < ^3(^2,0)}, 
[v 2 ,n 2 ],y < 0,Vi(x,y) < Vi(/x 2 ,0)}, 
[x ,v 2 ],y < 0,V 6 (x,y) < V 3 (u 2 ,y 3 )}, 
[vi,x Q ],y < 0,V 7 (x,y) < V^i,^)}, 
\Mi,vi],V <0,V 2 (x,y) <V 2 (m,0)}. 





= {x G 


n 2 


= {x G 


n 3 


= {x G 




= {x G 


n 5 


= {x G 


n 6 


= {x G 


n 7 


= {x G 


n 8 


= {x G 




vi 



Figure 2. 

Here yi > 0, y 2 > 0, y 3 < 0, j/4 < are solutions 
of the square equations 

V 1 (v 1 ,y 1 ) = V 1 ( l ii,0), 

V3(v2,y2) = F 3 (/x 2 ,0), 

Vi(v 2 ,y 3 ) = Vi(n2,0), 

V 2 (u 1 ,y 4 ) = V 2 (fii,0) 

For the derivatives of the functions Vj(x, y) along 
the solutions of system (24) we have the following 
relations 

X 

Vi = -2f(x)y 2 , V 2 = -2g(x) j f(z)dz, 

X 

V 3 = -2g{x) j f(z)dz, 

V'2 
X 

V 4 = -2g(x) j f(z)dz 

v\ 

X 

V 5 = ~2g(x) J f{z)dz - ey, 

V2 
X 

V 6 = -2g(x) J f(z)dz + ey, 

X 

V 7 = -2g(x) f f(z)dz + ey. 

v\ 

Therefore under the above assumptions concern- 
ing hi and (i 2 and for y ^ 0, x ^ Vj we have the 
inequalities 

Vi < on DiUfi 5 , V A < on Q 2 , 

V 5 < on fi 3 , V 3 < on fi 4 , 

l/ 6 < on fi 6 , V7 < on 7 , 

V 2 < on f2 8 . 



Note that by conditions 2), (20) for sufficiently 
small e we have inequalities 2/5 < 2/6 an d 2/7 > ys-i 
where 2/5 is a positive solution of the equation 

^4(^0,2/5) = V2(i/i,yi), 
2/6 is a positive solution of the equation 

^5(^0,2/6) = ^3(1/2,3/2), 
2/7 is a negative solution of the equation 

V 6 (x ,y 7 ) = V 3 (u 2 ,y 3 ), 
2/8 is a negative solution of the equation 

V 7 (x ,y 8 ) = V 2 (vi,y A ). 

(Fig. 2). 

Thus, we construct a transversal closed curve, 
which proves the first assertion of theorem. The 
instability of focal equilibrium gives the existence 
of cycle. □ 

Let us clarify the fact that the relations V4 < 
0, V5 < 0, V% < 0, Vr < are satisfied on the sets 
^3; ^6; ^7) respectively. 

Hold fixed the arbitrary e > 0, we choose \i\ an d 
fx 2 so much closer to a and to +00, respectively, 
that the minimal values of \y\ on the intersection 
of the closed curve (Fig. 2) and the band {x £ 
[^i, v 2 ]} are more than 



and 



- max 2 

e xe[u 1: i/ 2 ] 



— max 2 

e xe[vi,v2] 



X 

g(x) J f(z)di 



g( x ) / f(z)dz 



This implies the required inequalities Vj < 0. Thus, 
we have here [ij = IJ<j( £ ) an d 

lim fii(e) = a, lim fj, 2 (e) = +00. 

Represent now conditions 1) and 2) in terms of 
quadratic system (1). 
We have 

h 

Without loss of generality, we put 61 > 0. Condi- 
tions 1) and 2) are satisfied if 



< 2c 2 < 61, 
bi 

ai(2c 2 - 



A >0, 



(21) 



ai(6i&2 - aic 2 ) 
h\ >02 ' 

Besides, from conditions (21) it follows a posi- 
tive invariance of the half-plane {x > a} for qua- 
dratic system (1). 

Note that here c\ = and the parameters a 2 
and /?2 do not enter into conditions (21). 

A set of closed piecewise-smooth transversal curves 
can be constructed similarly. Therefore we have 
the following 

Theorem 4. Let conditions (21) be valid. Then 
any solution, of system (1) with initial data such 
that x(0) > a, tends as t — > +oo to a bounded 
attractor, placed in the half-plane {x > a}. 

Theorem 4 permits us to localize the search of 
the limit cycles of quadratic system (1). We re- 
mark that if in the half-plane {x > a} we have 
a unique unstable by Lyapunov focal equilibrium 
of system (1) and for it conditions (21) are valid, 
then system (1) has a periodic solution, placed in 
the half-plane {x > a}. 

Besides, under these assumptions we have f(x) > 
for x < a. Then, using the Lyapunov functions 

X 

V(x,y) = y 2 + J g{z)dz 

2a 

we can prove that a solution of system (1) with 
the initial data x(0) < a either tends as t — > +oo 
to the equilibrium, either to infinity, either leaves 
in a definite time the half-plane {x < a}. 

Using the above and Theorem 4 we 

Theorem 5. Suppose, conditions (21) are valid 
and in the half-plane {x > a} system (1) has a 
unique unstable by Lyapunov focal equilibrium. Then 
any solution of system (1) with initial data such 
that x(0) < a either tends as t — > +oo to the equi- 
librium, either to infinity, either to a bounded at- 
tractor, situated in the half-plane {x > a}. Any 
solution of system (1) with initial data such that 
x(0) > a tends as t — > +oo to a bounded attractor, 
situated in the half-plane {x > a}. This attractor 
has at least one cycle. 

Example 2. We put b± = f3\ = 1, a.\ = 0, a 2 = 
b 2 = -1, ci = 0, c 2 = 1/4, ai = -1, /3 2 = 2, a 2 = 
-1000. Here conditions (21) are fulfilled, g(x) ^ 
for x ^ and x > —1, equilibrium x = y = is 
unstable focus. Therefore system (1) has at least 
one cycle by theorem 5. 

The conditions of Theorem 5 effectively select 
in the space of parameters of system (1) a set, 
of positive Lebesque measure, in which the cycles 



exist. What the obtaining of the results of such 
kind is of present interest is remarked in [Arnol'd, 
2005]. 

V.A.Arnold (see [Arnol'd, 2005]) writes: "To es- 
timate the number of limit cycles of quadratic vec- 
tor fields on a plane, A.N.Kolmogorov distributed, 
as a mathematical practice, few hundreds of such 
fields (with the randomly chosen coefficients of 
quadratic expressions) among few hundreds stu- 
dents of the mechanics and mathematics faculty 
of Moscow State University. 

Each student had to find the number of limit 
cycles of his field. 

The result of this experiment was perfectly sud- 
den: all fields had none of the limit cycles! 

For a small change of coefficients of field, a limit 
cycle is preserved. Therefore in the space of co- 
efficients, the systems with one, two, three (and 
even, as was proved later, four) limit cycles make 
up open sets and if the coefficients of polynomi- 
als are drawn at random, then the probability of 
hitting in this sets are positive. 

The fact that this event did not occur suggests 
that the above-mentioned probabilities are, evi- 
dently, small." 

But the set of existence of cycles with selected 
by theorem 5 is not small. 
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